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Abstract 

We show that the index of BPS bound states of D4, D2 and DO branes in IIA 
theory compactified on a toric Calabi Yau are encoded in the combinatoric counting 
of restricted three dimensional partitions. Using the torus symmetry, we demonstrate 
that the Euler character of the moduli space of bound states localizes to the number 
of invariant configurations that can be obtained by gluing DO bound states in the C 3 
vertex along the D2 brane wrapped P 1 legs of the toric diagram. We obtain a geometric 
realization of these configurations as a crystal associated to the extra bound states of 
DO branes at the singular points of a single D4 brane wrapping a high degree equivari- 
ant surface that carries the total D4 charge. We reproduce some known examples of 
the partition function computed in the opposite regime where DO and D2 charge are 
dissolved into D4 flux, as well as significantly generalize these results. The crystal rep- 
resentation of the BPS bound states provides a direct realization of the OSV relation 
to the square of the topological string partition function, which in toric Calabi Yau is 
also described by a theory of three dimensional partitions. 



1 Introduction 

The recent conjecture of Ooguri, Strominger, and Vafa [T| that the indexed entropy 
of M = 2 black holes, in a mixed ensemble, factorizes, perturbatively to all orders for 
large charges, into the square of the topological string, together with subsequent checks 
for the case of IIA compactified on local Calabi Yau in [2] [3] jl] has generated much 
renewed interest in the structure of the partition functions of twisted M = 4 Yang-Mills 
living on D4 branes. In this work, we will discover a situation in which this theory 
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localizes to a discrete sum over equivariant configurations, associated to a new type of 
melting crystals, related to those of [5] which describe the topological A-model in the 
same geometry. Taking the limit of large D4 brane charge automatically reproduces 
the square of the A-model crystals, at the attractor values of the Kahler moduli. 

Consider D4 branes wrapping various 4-cycles of a toric Calabi-Yau, with chemical 
potentials for D2 and DO branes turned on. The partition function, studied in [0], is the 
generating function of the number (really index) of BPS bound states, which is given 
by the Euler characteristic of their moduli space. This can be computed in the A^ = 4 
worldvolume U{Ni) Yang Mills on the 4-cycles wrapped by iVj D4 branes, with some 
bifundamental couplings along the intersections. We will find that, at least locally in 
a single vertex geometry, the number of bound states is captured by the combinatoric 
counting of three dimensional partitions subject to a certain truncation. 

From the perspective of the remaining four dimensional spacetime, these bound 
states appear as BPS black holes, preserving half of the N = 2 supersymmetry, pos- 
sessing a classical horizon with nonzero area in the regime of large charges. The Kahler 
moduli live in vector multiplets in IIA theory, hence they are driven to their attrac- 
tor values at the horizon, which are determined entirely by the brane charges, via 
the quantum corrected prepotential. The Legendre transformation leading to a mixed 
ensemble, with chemical potentials (jp and </> a for the DO and D2 charges, results in 
attractor values that are independent of the prepotential, 

47T 2 i 

9to P = To : — n> ta = 9to P (p a + -<p a ), (1) 

as shown in jL. The BPS partition function in this ensemble factorizes in the large 
charge limit according to the OSV relation, 

z BP s(^,r;p ,P a ) = \Ao P (gto P ,t a )\ 2 , (2) 

to be interpreted as a perturbative expansion in the topological string coupling, g top ~ 
1/0°. We will be concerned with the case of vanishing D6 charge, p° = 0. 

The torus symmetry of a toric Calabi Yau induces an action on the moduli space 
of branes, enabling us to apply the equivariant localization theorem that the Euler 
character of the full moduli space is the same as that of the space of invariant configu- 
rations. Therefore we can assume that the DO branes wrap the fixed points, and the D2 
branes wrap fixed curves, which are the P 1 legs of the toric diagram. Imagine cutting 
the Calabi Yau along the toric legs, so that it looks like copies of the local C 3 vertex 
glued along P-'-'s over which the manifold is fibered as 0(—p) © 0{p — 2) for an integer p 
sometimes called the framing. Then the bound states of D2 branes can be recovered by 
matching equivalent configurations at the cut, and DO branes are localized to a single 
vertex. This is very reminiscent of the toric localization [7] of the topological string 
amplitudes, particularly as computed by the quantum foam D6 brane theory |Hj. 

The generic points in the moduli space of very ample D4 branes can be reached 
by turning on VEVs for the M = 4 adjoint field describing deformation in the normal 
direction inside the Calabi Yau, leading to a configuration of a single D4 brane wrapping 
a complicated smooth surface, C C X, whose homology is fixed by the brane charges. In 
this phase, the bound states that were described in terms of D2 and DO branes dissolved 
into flux of the Y\ U(Ni) quiver theory of intersecting D4 branes, become U(l) flux 
configurations associated to the large number of new divisors in the deformed surface. 
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The "D2" bound states wrapping curves in H2(C) > H^X) which are home-logically 
trivial in the ambient 6 dimensional geometry carry no D2 charge, and replace the DO 
flux bound states present in the Y\ U (iVj) phase. 

It is easy to see that such smooth D4 surfaces beak the (C*) 3 symmetry of the 
toric Calabi Yau. The crystal construction arises when we consider nilpotent VEVs 
for the adjoint field, which allows the use of the underlying toric symmetries to solve 
the theory. This corresponds to the limit in which the new divisors in C shrink to zero 
size, giving rise to the nontrivial DO bound states at the singular points that will be 
counted by the crystal. 

The bound states of D2 and DO branes to a single stack of N D4 branes wrapping 
a noncompact cycle in a toric Calabi Yau was studied in [3], where the theory was 
reduced to q-deformed Yang Mills in two dimensions. The D2 bound states wrapping 
a P 1 are parameterized by irreducible representations of U(N), which are associated 
to the holonomy basis of the Hilbert space of the two dimensional topological gauge 
theory. These configurations can be readily understood in the Higgs phase of the 
D4 theory, describing the generic point in the moduli space of normal deformations 
where there is a single D4 brane wrapping a complicated surface. In that case, the D4 
worldvolume is a N-sheeted fibration over the P , and hence, locally, the transverse 
motion of k bound D2 branes is described by k points in the fiber, with simply looks 
like N copies of C. Note that for the purpose of this discussion we can ignore the 
global obstruction to deforming a curve inside the D4 worldvolume. This is because 
the partition function separates into propagator and cap contributions, and the D2 
bound states can be equivalently analyzed in the propagator geometry, over which the 
bundles are trivial. Counting the number of ways of distributing k objects among N 
choices gives the relevant U(N) representations, associated to the moduli space 

Sym k (NC) 
Sn 

We will see this feature more clearly by choosing a toric invariant D4 worldvolume in 
which to perform the analysis, where these bound states will be described by certain 
torus equivariant ideals. 

Furthermore, there is an interesting and suggestive structure to the generating 
function of bound states in the vertex geometry of a stack of N D4 branes with only 
DO branes. In the q-deformed Yang Mills language, this is the vacuum cap amplitude, 
that is, the amplitude with trivial holonomy, which has been calculated in [3] to be 

^ w 4EHv- w =^ 4 n^ n y-v w 

' o-eSjv n>0 V q ' \<i<j<N 

The BPS states in which the DO charge is dissolved into U(N) flux are captured by 
the q-deformed Yang Mills amplitude, while the r] function counts one state for each 
fc-tuple of point-like DO branes bound to one of N identical D4 branes, noting that the 
effective Euler characteristic of the cap is x = 1- This amplitude can be expressed as 

n>0 v H ' j=l j=l y H 1 j>N y y ; 

associated to bound states of fc-tuples of DO branes with k ground states, truncated 
at N. This function nicely interpolates between the 1/rj for a single D4 brane, and 
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asymptotes to the McMahon function, which counts bound states to a D6 brane. The 
truncation has a transparent interpretation in terms of the crystal we will now proceed 
to describe. 

Consider the generic situation of N, M, and K D4 branes wrapping the three in- 
tersecting 4-cycles in the local C 3 vertex, with bound D2 branes in the legs, and n DO 
branes at the vertex. We will first find that the equivariant BPS bound states of p D2 
branes wrapping the intersection curve of N and M D4 branes are partitions of p that 
do not contain the point (M + 1, N + 1). In the special case M = 0, this reduces to 
Young tableaux with ./V rows, which are related to the U(N) representations that span 
the q-deformed Yang Mills Hilbert space. Such bound states have D2 charge p, and 
induce some DO charge which depends on the local bundle over the wrapped curve via 
the intersection form. Fixing such Young diagrams, R, Q, and P, as the asymptotic 
D2 brane configuration in the vertex geometry, we will show that each contribution 
to the index of BPS states carrying an additional n units of DO brane charge can be 
associated to a three dimensional partition of n in the truncated region of the octant 
obtained by deleting all the points beyond (iV + 1, M + 1, K + 1), with asymptotic 
behavior R,Q,P, as shown in figure 4. 

The organization of this paper is as follows: In section 2, we review the theory of 
BPS D4 brane bound states, particularly in the local toric geometries studied in |2j [3] 
and and the OSV relation to topological strings pQ. In section 3, we will examine 
the worldvolume theory on these branes in the Higgs branch. Giving a nilpotent VEV 
to the adjoint field, which controls normal deformations inside the Calabi Yau, reduces 
the partition function to a sum over toric invariant ideal sheaves of certain singular 
surfaces. In section 4, we study the bound states with D2 branes in the toric "leg". 
In section 5, we compute the index of DO bound states in the vertex geometry. In 
section 6, we see that the OSV factorization is manifestly verified in all cases where 
these methods apply. Finally, in section 7, we make some concluding remarks. 

2 BPS bound states of the D4/D2/D0 system 

The remarkable conjecture of OSV relates the large charge asymptotics of the in- 
dexed degeneracy of BPS bound states making up supersymmetric black holes in four 
dimensions with the Wigner function associated to the topological string wavefunction. 
Consider compactifying IIA theory on a Calabi Yau manifold, X, which gives rise to an 
effective TV = 2 supergravity in the remaining four dimensions, with vector multiplets 
associated to the Kahler moduli. The microstates of \ BPS black holes in this theory 
are bound states of D6/D4/D2/D0 branes wrapping holomorphic cycles in X. 

The effective theory in four dimensions has U(l) gauge fields obtained by integrating 
the RR 3-form, C 3 , over 2-cycles, [D a ] G H 2 (X) for a = 1, . . . h l ^(X), in the Calabi 
Yau. The D2 branes wrapping these 2-cycles and D4 branes on the dual 4-cycles, 
[D a ], carry electric and magnetic charges respectively under these C/(l)'s. There is one 
further U(l) vector field, arising from the RR 1-form in ten dimensions, under which 
the D0/D6 branes are electrically (magnetically) charged. 

We will investigate this partition function at D6 charge in the local context of 
toric Calabi Yau, where the theory on the branes is given by twisted AT = 4 Yang 
Mills on each stack of D4 branes, coupled along their intersections. The action of the 
U(N) worldvolume theory of N D4 branes is BRST exact, and turning on the chemical 
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potential for bound D2 and DO branes inserts the closed (but not exact) operators, 

into the action. These topological terms count the induced D2 and DO charge of a 
given gauge connection. 

This Vafa-Witten theory jS] computes the indexed degeneracy of BPS bound states 
with partition function equal to 

Z YM (Q 4 a ,<p a ,<p°)= Y, n(Qt,Q 2a ,Q ) exp {-Q <p° - Q 2a ^ a ) , (7) 

Q2 a,Qo 

in the mixed ensemble of fixed D4 (and vanishing D6) charges. 

A collection of D4 branes in the OSV mixed ensemble will be well described by 
semiclassical ten dimensional supergravity as a black hole with large horizon area when 
the branes live in a homology class deep inside the Kahler cone Given D4 branes 
wrapping such a very ample divisor, [D], one can show that the attractor values of the 
moduli are large and positive, 

me(t a ) » 0, (8) 
as well as the leading classical contribution to the black hole entropy, 

C abc t a t b t c > 0, (9) 

since the intersections numbers of [D] with 2-cycles are positive. It is in this regime 
where we expect the brane partition function to factorize into the square of the topo- 
logical string wavefunction. Moreover, these are precisely the conditions on the 4-cycle, 
with homology [D] that is wrapped by the D4 branes, to have holomorphic deforma- 
tions. The existence of such global sections of the (2, 0) bundle over the 4-cycle in 
the compact case would allow us to turn on a topologically exact mass deformation 
and solve the theory following |3J. For branes wrapping noncompact 4-cycles, the 
mass deformation results in the reduction to a theory on the two dimensional compact 
divisors. 

This two dimensional theory was shown to be the q-deformed U(N) Yang Mills 
on the divisor of genus g in [5], who then solved the theory using techniques of 
topological gauge theory in two dimensions, obtaining the result 

Z YM = ^(dim 9 ^) 2 - 2 ^I^W e ie|^| 5 (1Q) 
n 

where 1Z is an irreducible U(N) representation. The more general situation of D4 
branes intersecting in the fiber over the intersection point of two curves was solved in 
the toric context by [I] , where the q-deformed Yang Mills theories found on each curve 
were coupled by operators inserted at the intersection point, obtained from integrating 
out bifundamentals. In this work, we will check the results of the crystal partition 
function against these examples, as well as generalizing to configurations of D4 branes 
intersecting along the base. 

The BPS partition function in the large charge limit was shown to factorize, con- 
firming the generalization of the structure Z ~ |^to p | 2 to the noncompact case with an 
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extra sum over chiral blocks parameterized by SU(oo) Young tableaux, P, P' control- 
ling the noncompact moduli. In particular, for the case of N D4 branes in the Calabi 
Yau 0(-p) 0(p - 2) -> P 1 studied in the large N 't Hooft limit gives 

Z = E E Zf^ + (t+pg s i)Zf^-(i-pg s ll (11) 
t& P,P> 

where the topological string amplitudes are 

xY j q E ^ KR e- tlRl W PR (q)W P n R (q) (1 2 ) 
ZfP'-{t) = {-ir\ + \ P '\zf t f,t{t), 

and 

4vr 2 p-24tt 2 <b 2 

*-*T ( = — ^" 2 'V (13) 

are the attractor values of the moduli. 

Recall that in toric Calabi Yau all compact 4-cycles have positive curvature, and 
thus negative self- intersection number. Branes wrapping these rigid cycles will con- 
tribute negative real part to the attractor moduli. It is possible, however, that even 
in this rigid case a modified local analysis is still possible, in which an additional con- 
straint for each wrapped compact 4-cycle is introduced on the sum of local bundles over 
its toric divisors. For example the twisted = 4 Yang Mills on the rigid divisor P 2 in 
the Calabi Yau 0{— 3) — > P 2 is solved in [H], ^U], and These partition functions 

do not factorize, and correspond to the regime in which nonperturbative baby universe 
effects are not suppressed ^2] (!]. 

More applicably in the context of OSV, some of the branes might wrap rigid sur- 
faces, as long as the total charge is the in Kahler cone. In this case one can still expect 
a factorization into the topological and anti-topological string, and a simplification of 
the Yang Mills theory on the compact rigid surfaces may occur due to the intersection 
with sufficiently large numbers of branes wrapping very ample divisors. 

Another interesting feature of the D4 brane partition functions is their behavior 
under s-duality of the twisted theory in four dimensions. This is not manifest in 
the crystal formalism, since in our situation of noncompact D4 branes, the boundary 
conditions at infinity in the noncompact directions transform in a complicated way. It 
would be interesting to understand more completely how the duality is realized in our 
picture, although one aspect is already clear. Namely, we are really counting the vector 
bundles associated to the dual flux of the topological Yang- Mills theory. This explains 

4tt 2 

why our expansions appear most naturally in terms of q = e~ 9s = e ^° rather then 
q = . We will also use the variable 9 a = ( / ,a - 

The boundary conditions that are naturally associated to the crystal description 
are that no bound D2 branes should reach oo in the noncompact D4 worldvolume. This 
turns out to be the same condition used in the derivation of q-deformed Yang Mills [2] 
PI) but is the opposite of the standard "electric" convention from the four dimensional 
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worldvolume point of view, which would require the holonomies of the gauge field to 
vanish asymptotically. These are exchanged by the action of s-duality. For exam- 
ple, the instantons contributing to the U(l) Yang Mills amplitudes on 0{— p) — > P , 
expressed in terms of the gauge fields, are given by the lift of line bundles on the 
P , with noncompact toric divisors lying in the fiber over the fixed points. The dual 
configuration has the P 1 itself as toric divisor, and it is this that we shall count. 

2.1 Turning on a mass deformation for the adjoints 

The deformations of the wrapped surface are encoding in the adjoint fields of the 
twisted Af = 4 U(l) gauge theory living on the D4 brane. The torus symmetry 
of the underlying geometry can be used to localize the D4 brane to the equivariant 
configuration, by the usual method of turning on a Q-trivial mass deformation, 

W = mUV + wT 2 , (14) 

parameterized by the choice of u, a global section of the (2, 0) bundle over C, where 
U and T are scalar adjoint fields associated to motion in spacetime, while V governs 
deformation normal to C inside the Calabi-Yau. The superpotential, W, must live in 
a (2, 0) bundle over C. 

In our case, there is an important new feature to this story, due to the compactifi- 
cation of the normal direction in the Calabi-Yau. The adjoint field, V, will take values 
in some compact space, which, in general, would greatly complicate the analysis. As 
explained in t 4 this issue can be avoided when the D4 branes wrap noncompact cycles, 
since one is free to choose boundary conditions at infinity corresponding to a partic- 
ular surface in the given homology class. Changes of these boundary conditions are 
determined by non-normalizable modes in the four dimensional theory on the branes, 
and we should not sum over them, at least in the local case. It is sensible to choose 
this divisor to be torus invariant, for consistency with the toric symmetry of the local 
geometry (in fact, it is only in the equivariant sense that we can even appropriately dis- 
tinguish various 4-cycles which would otherwise appear homologous in the noncompact 
Calabi-Yau). 

It is gratifying to understand this point in greater detail, for concreteness in the 
example of local P 1 x P. Let us consider a single D4 brane wrapped on the 4-cycle, 
C, given by 0{— 2) — > P , whose normal direction is a compact P . This seems to 
imply that the field, V, should be P 1 valued, however this is modified by the fact that 
our 4-cycle, C, is nontrivially fibred over this P 1 . The local coordinate normal to the 
brane is thus a section of the 0(—2) bundle, rather then the trivial bundle. Globally, 
this implies that equation (|14|) will have two saddle points where this section vanishes. 
In order to preserve the torus symmetry, we must choose an equivariant section of 
the 0{—2) bundle, which implies that the mass deformation localizes the D4 brane to 
exactly the two invariant surfaces with the homology class of [C], namely those sitting 
above the north and south pole of the other P , exactly as expected. 
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3 Localization in the Higgs branch of the theory 
of intersecting branes 

The worldvolume theory of D4 branes can be localized, due to the topological nature 
of the partition function, to a computation of the Euler characteristic of the moduli 
space of anti-self dual instanton solutions. This is because the action of the twisted 
theory is exact in the cohomology of the twisted supercharge, Q, and the chemical 
potentials for DO and D2 branes couple to the instanton charge, J F A F, and Jj^j F 
respectively, which are topological invariants. It can be shown that these anti-self dual 
connections exist uniquely for each stable U(Ni) vector bundle (or certain sheaves in 
a more general, singular context). 

It is also necessary to take into account the contribution of additional bound DO 
branes, which are given by sheaves with singularities at points where the DO branes 
are located. This gives rise to a universal contribution of r/~ x of the eta function 
to the power of minus the Euler character The eta function appears as the 

generating function of Young tableaux associated to the Hilbert scheme of points in 
two dimensions. Such effects will be automatically included in our case, where all of 
the bound states will have to be described exclusively in the sheaf language, since the 
D4 branes will wrap singular surfaces. 

Any collection of intersecting D4 branes whose total homology class is very ample 
can be deformed in the normal directions, generically producing a single D4 brane 
wrapping a complicated surface. Thus the quiver of U(Ni) twisted Af = 4 Yang Mills 
theories with bifundamental interactions can be replaced by the £7(1) gauge theory on 
this high degree 2-fold, C, which computes the same partition function. In both cases 
the adjoint field associated to motion of the brane in the Calabi-Yau directions will 
generally be compactified in a highly nontrivial manner. 

For example, consider wrapping N D4 branes on the cycle 0{— 3) — ► P 1 in the Calabi 
Yau 0(-3) 0(1) -> P 1 , and M on the rigid 0(1) -> P 1 cycle. Then the attractor 
value of the Kahler modulus, 9\t ipi = 4f (TV — 3M), will be positive for N > 3M, hence 
we should be able to write down a smooth surface with these charges. The projective 
coordinates are given by (z,v;x,y) ~ (Xz, Xv; Xx, X~ s y), where z = v = is deleted. 
The deformed 4-cycle is described an equation of the form 

x N y M + ax N ^y M - 1 + ■■■ + bz N ~ 3M + ■■■+ cv N ~ 3M = 0, (15) 

where each term has degree N — 3M > 0. 

Therefore we really want to calculate the generating function for the Euler char- 
acters of the moduli space of ideal sheaves (ie. rank 1 stable sheaves) on C. To do 
this, we will extensively use the toric geometry of the underlying Calabi Yau manifold. 
For such techniques to be applicable, the surface, C, must be chosen carefully to be 
invariant itself. For a high degree surface with positive self-intersection of this type, 
the only possibility is a highly singular realization as a "thick" subscheme living on the 
original intersecting toric divisors. This is achieved in the worldvolume gauge theory 
by giving a nilpotent VEV to the adjoint field V that controls the normal deformations 
inside the Calabi Yau. 

The ideal sheaves on such an affine space are equivalent to ideals of the algebra 
of functions, and can be counted in the equivariant setting as described below. By 
cutting the subscheme C into pieces, we can obtain the full answer by gluing vertices 
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via propagators, analogously to 8 . Note that toric localization is crucial for this 
procedure to work, with DO branes localizing to the vertices, and D2 branes to the 
legs. The vertex amplitudes, associated to bound states of DO branes with the D4 
and D2 branes, in the chemical potential, will turn out to be nicely encoded as the 
statistical partition function of a cubic crystal in a restricted domain. 



3.1 Solving the twisted J\f = 4 Yang Mills theory 

The key point is that we have replaced the complicated system of intersecting stacks of 
Ni D4 branes described by bifundamental-coupled U (iVj) Yang Mills theories by a single 
D4 brane wrapping a nonreduced subscheme of X. The U{\) theory on that brane can 
be solved using toric localization, and the well known fact that the Vafa-Witten theory 
computes the generating function of stable bundles. The partition function is given by 

Zym= £ q n e^ e « A «x(M n , [A] ), (16) 

n£H°,[A]eH 2 

where .M,^^ is the moduli space of coherent sheaves, n is the instanton number (equiv- 
alently, the DO charge), [A] is the D2 charge, and a indexes H.2- The effective number 
of bound states with given charges is given by the Euler characteristic of the instanton 
moduli space. The brane charges are given exactly by the coefficients of the exponential 
Chern character, which is, in our situation, 

1 , 

q D0 = 2ch 3 = c 3 - C!C 2 + -Cj 
1 2 

QD2 = ch 2 = -c 2 + ~Cj ( 17 ) 
q m = ch] = [C] 
QD6 = 0, 

where we are ignoring the gravitation contribution from the first Pontryagin class of 
the Calabi Yau, which gives an overall contribution independent of the D2/D0 bound 
state, and is in fact somewhat ambiguous in the noncompact setting. 

We want to determine the contribution of DO branes bound to the vertex in the 
toric geometry the looks locally like C 3 . In order to be able to apply the machinery of 
equivariant localization, we must choose a torus invariant representative of the homol- 
ogy class \C\. This is impossible for smooth surfaces, thus we are forced to consider the 
singular, nonreduced subscheme described by the polynomial equation x N y M z K = 0, 
where N, M, K are the number of D4 branes wrapping the equivariant homology classes 
associated to the three planes in C 3 . Because the worldvolume is so singular, we need 
to be careful when determining the moduli space of line bundles. Looking at the ring 
of functions on the surface, it is clear that it will approach the algebra 



_ C[x,y,z] 
(./•V':'- 



as we deform a smooth surface to the singular limit. This corresponds to a non-reduced 
subscheme of C 3 , since the polynomial x N y M z K does not have separated roots, which 
translates the fact that the adjoint fields of the gauge theory have nilpotent, non- 
diagonalizable VEVs. 
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Given any very ample line bundle on the affine scheme C, we can associate to it in 
the usual way an ideal of A by considering its sections. For a general line bundle, £ , 
corresponding to asymptotic representations with negative row lengths along a given 
curve, we first take the tensor product with a line bundle, C n , over the D4 worldvolume 
of sufficiently positive first Chern class, n, to obtain a very ample bundle. The bundle 
C n must descend from a bundle, 0(n), on the physical 4-cycle, namely, its divisor is 
simply the multiply wrapped curve in question. The instanton number can easily be 
calculated to be 



since ch{£) <g> ch(C) = ch(£)ch(C). Hence we can effectively concentrate on the SU(N) 
component of the partition function, from which the full answer can be trivially recon- 
structed by summing over the choice of such line bundles, C n . 

Now we can use the action of the toric symmetry to considerably simplify the 
problem. There is a equivariant localization theorem which relates the cohomologies 
of the moduli space A4 n \m of bundles (or rather sheaves) on C with Chern classes 
n E Hq = 7L and [A] G H2 to those of invariant sheaves, 



that is, the Euler character of the moduli space is captured by the considering only the 
equivariant bundles. In more physical terms, the partition function involves summing 
over BPS states of D2 branes wrapping the legs of the toric diagram, that is the torus 
invariant P^s, and DO branes at the equivariant points, which are the vertices, bound 
to the D4 brane wrapping the nonreduced cycle C. 

In the next section, we will determine the number of equivariant bound states with 
wrapped D2 branes, and compute the induced DO charge. They will correspond to a 
choice of representation along each leg, with rank determined by the D4 charges (for 
example, we will have U(N) representations for D2 branes bound to a stack of N D4 
branes). By cutting the Calabi Yau along the legs of the toric diagram, to effectively 
obtain glued C 3 geometries. Clearly this cutting of the Calabi Yau is only valid in 
the toric context, since we have seen that the sum over D2 branes localizes to those 
wrapping only the legs. The situation here is analogous to that studied in [7j where 
the topological string amplitudes on toric Calabi Yau were shown to reduce to gluing 
copies of the C 3 vertex, because the worldsheet instantons localized to equivariant 
curves. Similar localization arguments were used by (Hj in the more precisely related 
context of counting instantons of the U(l) gauge theory living on a D6 brane, and 
we will not repeat them here. This defines a vertex, Vrqp, depending on up to three 
intersecting stacks of D4 branes intersecting in C 3 , and the asymptotic representations, 
R,Q,P, associated to the D2 configuration. Gluing these vertices by the propagators 
along the legs, we obtain the full partition function. 

In the vertex geometry, the torus invariant bound states of DO branes to the D4 
worldvolume in the nilpotent Higgs phase are given by the equivariant ideals of the 
algebra (|18|) . The toric symmetry acts by (C*) 3 multiplication on the coordinates, 
(x, y, z) — ► (Aix, A22/, A32;), hence the only invariant ideals are those generated by mono- 
mials. Depicting the monomials, x n y m z k , of the polynomial algebra, C[x,y,z], as the 
lattice of points, (n,m,k), in the positive octant, the algebra relation x N y M z K = is 



ch 2 (£ ® C n ) = ch 2 {£) + ci(C n )ci(£) + -cj(£ n ) 



(19) 




(20) 
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represented by deleting all points (monomials) set to zero. The remaining points give 
a three dimensional partition, which has been "melted" from the corner of the cubic 
lattice. 

Consider a DO bound state described by a toric ideal, I<A, spanned by monomials, 
x n y m z k , for (n,m,k) E S. Clearly if x n y m z k £ X, then so is x n+1 y m z k and so on, 
hence the quotient algebra ^ wm be depicted as a three dimensional partition, in the 
truncated region described by the D4 worldvolume algebra, A. This quotient is the 
natural algebra associated to the DO worldvolume, the codimension 2 analog of the 
divisor of the sheaf on C. Hence the DO charge of this bound state is exactly the 
dimension of the quotient algebra, which is the number of points in the partition. The 
region in which the configurations are restricted to live is defined by the polynomial 
relations of the algebra given by (|18|). shown graphically in figure 4. 



4 Propagators and gluing rules 

Typically, our D4 branes will wrap 4-cycles containing compact toric divisors joining 
two fixed point vertices. These geometries will always look locally like O(-p) — » P , 
for some p, possibly intersecting other D4 branes along the toric divisors in the Calabi 
Yau, X. For simplicity, we will first analyze this geometry in isolation, with a single 
stack of N wrapped D4 branes. They can be described in the above language as the 
subscheme, C C X, associated to the graded algebra 

A= C[X: J> V > Z] , (21) 

with weights defined by 

(x : t; y, z) ~ (Xx : At; X~ p y, X p ~ 2 z), (22) 

for A G C x . 

Rank N vector bundles on the 4-cycle can now be described as rank 1 sheaves 
over C, and thus by ideals of A. Moreover, invariant bundles whose corresponding 
equivariant divisors are (possibly non-reduced) multi- wrappings of the IP 1 are equivalent 
to homogeneous ideals of the form 

l=({y R >z i - 1 }), (23) 

for i = 1, ... TV. These ideals are in one to one correspondence with irreducible rep- 
resentations of U(N) with positive row lengths. More general representations can be 
obtained as explained above in equation (|19|). We proceed to calculate the Chern char- 
acters of such bundles, which will serve as the background configuration of our crystals. 
This might seem difficult, since the base manifold is now non-reduced, so we will need 
to relate these sheaves to ones living over the full Calabi- Yau geometry. 

In a similar manner, the propagator along legs where stacks of D4 branes intersect 
can be determined by counting equivariant ideals of the algebra 



C[x : t; y, z] 
= (y M z N ^ ' 



They will be of the form 

I=({y Ri z i - 1 }), (25) 
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where the Young tableau R lies in the region restricted by N and M, that is, < Ri < 
M for i > N. We will refer to such Young diagrams as being of type (N,M). The 
following analysis of the induced DO charge for rank N bundles is easily extended to 
include this case as well. 

The subscheme, C, is embedded via a map 

i:C^X, (26) 

allowing us to push forward a sheaf, £, to i*£ , a sheaf on X with support only along 
the 4-cycle. This corresponds to the fact that the ideals of A defined above are clearly 
also ideals of the structure algebra of X itself. The Chern characters of i*£ can be 
determined using equivariant techniques to be jSJ 

2ch s = \\R l \\ 2 + ^k e 

" " 2 (27) 

ch 2 = \R\iF 1 }. 

This is not quite the whole answer, since those Chern characters have contributions 
not only from the bundle, £, but from the embedding of the subscheme C itself. There- 
fore we must essentially subtract the corrections arising from the normal bundle to C. 
More precisely, the Chern character and the push forward do not commute, in a way 
measured by the Grothendieck-Riemann-Roch formula |14j . 

ch{U{£))=U{ch{£)-Td{T f )), (28) 

where / : Y — ► X, ch is the Chern character, Td is the Todd class, and Tf is the 
relative tangent bundle. Our situation of an immersion into a Calabi-Yau manifold is 
a significant simplification of the special case of the GRR formula considered in |15j . 
We find that the Chern characters are related by 

ch s (uS) = i* (ch 2 (£)) - i[C\ ■ i, (cfcj(f)) = ch 2 (£) - N^-\R\ 

ch 2 (i*£) = ncht{£) = \R\ [F 1 ], 

up to trivial terms that are independent of R. Applying this in reverse to the Chern 
characters from the three dimensional point of view written in (|27|). we conclude that 
the brane charges are 

Qdo = |«J2 + H^ll 2 + ^N\R\ = ^C 2 (R) + WR'W 2 - N\R\ 

QD2 = \R\- 

It will be useful later to separate these terms into a propagator depending on the 
bundle, O(-p), and a piece that will be included in the vertex amplitude calculated 
below. Therefore we find the toric propagator to be 

<ficm e ie\R\^ (31 ) 

as expected from q-deformed Yang Mills, and the remaining pieces of the second Chern 
character give the corresponding factors, 

exp{- 5s Q||^|| 2 -^|i?|)|, (32) 
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in the two vertices being glued. 

These bound D2 branes also contribute to the induced DO charge when they in- 
tersect at the vertices. In particular, the intersection numbers give terms of the form 

N 

Y, R iQi> ( 33 ) 

i=l 

for D2 branes bound to the same stack of iV D4 branes along two intersecting legs of 
the toric diagram. We will see in the next section that these terms exactly cancel a 
correction to the vertex associated to the number of boxes deleted by the background 
D2 branes, as expected. 

If there are two intersecting stacks of D4 branes, the induced DO charge can be 
similarly computed as 

Qdo = |||«H 2 + ^ll^ll 2 + P ~^N\R\ - P -M\R\, (34) 
where the D2 bound state is of the type found in (|25[). 



5 The crystal partition function 

We will use the transfer matrix method, described, for example, in [T^], to evaluate the 
partition functions of these crystal ensembles. Recall that cutting a three dimensional 
partition along diagonal planes, y = x + t, gives a sequence of interlacing Young 
diagrams, {v(t)}, indexed by t. In particular, 

u(t + 1) >- u(t),t < 0, (35) 

and 

v{t + 1) -< i/(t),t > 0, (36) 
where two Young diagrams, and u, are said to interlace, denoted by (jl ~< v, if 

V\ > 111 >V2>^2> ... (37) 

It is now convenient to use the correspondence between Young diagrams and states of 
the NS sector a free complex fermion (or boson, via bosonization) in two dimensions, 
to express the crystal partition function as a sequence of operators acting between the 
states associated to the asymptotic Young diagrams. 

Given a Young diagram, v, define its Frobenius coordinates as 

1,1 
a n = v n -n+-, b n = u n -n+-, (38) 

where the index runs from 1 to d, the length of the diagonal of v. Then the fermionic 
state, 

d 

l^>=IIOd0>, (39) 

n=l 

is nontrivial and uniquely associated to since the a n and b n are distinct and uniquely 
determine v. 
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The key point in this transfer matrix approach is that the operators 



r ± (z)=exp(£^^V (40) 

\n>0 / 

constructed out of the modes, J n , of the fermion current, V , *V'> generate the interlacing 
condition. That is, 

r+(l)|i/ >= >, and (41) 

r_(l)|i/ >=£>>. (42) 

The number of boxes in a single slicing is easily computed using the Virasoro zero 
mode, Lo, to be 

q L °\u >= q \"\\ v >, (43) 

which can be applied sequentially to evaluate the crystal action of a given partition. 

It is important to take into account a finite boundary effect arising from the differ- 
ence between the number of boxes falling inside a region bounded by diagonal slicings 
as opposed to a rectangular box. If we regulate the crystal partition function by plac- 
ing it in a rectangular box of length L along some axis, then for L sufficiently large, 
only the background asymptotic representation, R, will contribute near the boundary. 
We will count the total number of boxes in the three dimensional partition, and hence 
must subtract off the background contribution, q L \ R \. The diagonal slicing includes an 
additional number of boxes, which are easily counted to be 



2" 11 2 1 1 



1^.(^-1), (44) 



or its transpose, depending on the orientation of R as in [3]. 

This is still not quite right when there are asymptotic representations along more 
then one leg, since the subtraction of +-L2IQI +-£ , 3|-P| over counts the number of 
boxes in the background by the intersection at the origin. 2 For example, the intersection 
of two partitions with the parallel orientation has 

^RiQi (45) 

i 

points. It will turn out that this effect is cancelled by the contribution of the Chern 
classes of the background divisors themselves. 

Define the normalized vertex in terms of the crystal partition function, P(R a , N a ), 
to be, schematically, 



= exp 



J y l^2N a \R a \-^2\\R a \\ 2 + 2^2R a -R b -4R 1 ■ R 2 ■ R 2 \ 1 P(R a , N a ), 

I V a a a <b / ) 

(46) 

where a ranges over the nontrivial asymptotics, R a • R b is given by (|45[). the final term 
is the triple intersection, and transposes must be used appropriately depending of the 



2 This mistake has actually appeared frequently in the literature. 
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Figure 1: Shown in red is a diagonal slice of the crystal describing 5 coincident D4 branes 
in the vertex. The meaning of the representation v will be discussed in section 6. 



orientation of the Young tableaux. These terms are perfectly cancelled by the Chern 
classes of the background bundle, which we determined in (1321) . 

We will make extensive use of various identities involving the operators, T, and 
skew Scaur functions, which are collected in Appendix A. 

5.1 One stack of N D4 branes 

Let us apply these methods first to the simplest case of a single stack of D4 branes. 
Our analysis in this subsection will parallel that of Okuda |17j . where the special 
case of a single nontrivial asymptotic representation was discovered in the context of 
Chern-Simons theory. 3 

Regulating the crystal by placing it in a box of size Lx N x oo, we apply the transfer 
matrix method as shown in figure 1, in a background determined by the representation 
Q along the vertical leg, with asymptotic R to the left and • required by the truncation 
on the right. The pattern of interlacing operators is altered by the presence of Q, as 
explained in 0, thus we see that the normalized partition function defined by (fl6|) is 

S RQ = q %W+\Q\) q l\\QTHW q -I<W <R\f[ {q L °r ± (l)) q L ° ]] (r T (l)g L ») |0 >, 

n>0 m=l 

(47) 

where the pattern of pluses are minuses is determined by the shape of Q. Commuting 
the q Lo, s to the outside, and splitting the middle one in half, we obtain 

N 

^IIQ'II 2 <R\]J T±(q- n+1 2) Y[ r T (g m -!)|0 >, (48) 

n>0 m=l 

in the L — > oo limit. Next, we commute all of the r + to the right, where they act 
trivially on |0 >, noting that the commutator, Z, depends only on the representation 

3 In fact, the original motivation of this work was to understand the underlying structure of that relation! 
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Figure 2: Here we rotate the crystal with trivial asymptotics associated to the U(2) theory 
with no bound D2 branes. 



Q. This gives 

N 

^IIQ'II 2 <R\Y[ r_(g- Qm+m -^)|0 > Z(Q) = q^ Qt W 2 s R {q- Qm+m ^)Z{Q). (49) 

m=l 

To determine Z(Q), note that 

S. Q = glllOT-£l<JIZ(Q) = S Q . = Z{-)s Q {q-P"). (50) 
Therefore, pulling everything together, we find that 

S RQ =S.. s Q (q- pN )s R (q' pN - Q ), (51) 

which is exactly the Chern Simons link invariant obtained from q-deformed Yang Mills 
in and gj! 

For future reference, we will derive another useful formula for S R q by slicing the 
crystal with a different orientation as depicted in figure 2. Applying manipulations 
similar to the above, we find that 

S RQ = gK-^l^-^IOI+ll^l'+IIQII^g-^ll^ll'+IIQ'H 2 ) 

x < a*| H r+iq-^vj, n t^-^iq* >, 

n>0 m>0 

(52) 

where Vn {P%) ls the projection operator onto the space of representations with at 
most N rows (columns). This implements in the transfer matrix language the fact that 
the three dimensional partition has height truncated by N, by enforcing the constraint 
on the largest diagonal tableaux. Writing this in terms of Schur functions, we find that 

S RQ =S^-)\ R ^%-^ C ^+ C ^ £ s P/R { q O)s P/Q {q% (53) 

PeU(N) 

where we have used (f79*|) to remove the transposes. 
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ure 3: The crystal for iV = 2 intersecting M — 1 D4 branes is oriented as shown, where 
green planes are the boundaries of the allowed region. 



5.2 N D4 branes intersecting M D4 branes 

Suppose we have N D4 branes wrapping the x — y plane, M in the y — z plane, with 
asymptotic conditions defined by a U(N) representation, R, along the x-axis, and 
a U(M) representation, Q, in the z direction. Then we will proceed to check that 
the D4 brane partition function is given by the vertex Vrq of the statical ensemble 
of three dimensional partitions restricted to lie within the volume of the D4 branes 
with asymptotics R and Q, normalized according to (|46|) as before. Using the transfer 
matrix method, with diagonal slicings oriented in the planes y = x + n, indexed by n 
(see figure 3), we obtain 



V (N,M) ||Qt||2_^|ij|_f |q| L |i?| -l||i?|| 2 +i|i?| x 

R,Q " " " 

L-M M ( 54) 

< R l \ II (9 L °r+(i)) q Lo K [J (r±(i)^) n (v^q^) |o >, 

n=l m=l fc>0 

where Q determines the pattern of T±, as before, and the position of the projection, 
Vjy, is fixed by the truncation constraint that the point (N + 1, M + 1, 1) is not in the 
partition. We commute the Lo's to the outside, dividing them at the plane y = M + z 
where we enforce the projection, to obtain in the infinite L limit, 



-M|i*|-ic 2 (fl)+i||Q*|| 2 -f \Q\ 



2-^-^2,,^ I, — 2-|v*| ^ < R l \ Y\ T + {q~ n+ ^)\A l >x 
AeU(N) n>0 

M 

< a i \ n r±(<r-5) n r T(9 M+fc -^io > . 

m=l fc>0 

(55) 

Now, commuting the remaining r + to act trivially on |0 >, we pick up a factor, 
Y(Q), independent of R. Writing the resulting expression is terms of Schur functions, 
we have shown that 

V«Q = g- M ' i? l-^ 2 ^ + ^l Qt H 2 - i 2 £ ' Q 'T(g) £ s At/Rt ( q -0)s At (q- Qt - p )« MlAl - (56) 

AeU(N) 



17 



Note that the commutator term, Y(Q), can be expressed as 

M 

<o\l[r ± ( q ™-h)l[r T (q M + k - 1 z)\o>, (57) 

m=l k>0 

which, up to a shift in the origin of slicings by M units, is exactly the M-truncated 
crystal with asymptotics Q and •. Thus the results of the previous subsection imply 
that 

Y(Q) =qf ISIg-HW SQ ( g PM). ( 58 ) 

Manipulating our expression with the Schur function identities ()79|) and (|8U|) . we 
find 



/•« \ AeU(N) j 



(59) 

x {s Q (q^)^(q Q+PM )) ■ 
We can now recognize the terms in parenthesis in (|59|) to be exactly the Chern- 
Simons link invariants, g^W+^G*) sgj (g) and S^fa -1 ) = S^. Therefore, re- 
calling that C^ N) (P) = C^ M) (P) + (AT - M)|P| for a J7(M) representation P, we find 
that 

V*, = E ^M^H^ 1 ^, (60) 

Pel/(M) 

noting that Spg^ vanishes unless P is a U(M) representation. This is in complete 
agreement, up to overall U(l) factors that were not carefully taken into account, with 
the result of [1], where this partition function was derived by reducing the D4 brane 
world volume theories to q-deformed Yang-Mills in two dimensions, coupled by bifun- 
damental fields living at the intersection point. 

Thus far we have only considered the case of trivial asymptotics along the intersec- 
tion plane of the branes, which is natural assuming that direction is noncompact in the 
toric Calabi-Yau. Lifting that restriction, we proceed much as before after replacing 
|0 > by \P l > in ([HK]). to obtain 

-M|ii|-i||Rl| 2 -i||P[|2+f \Q\-±\\Q*\\*(\\R\ 



1 X 



w \AeU(N) J 

(61) 

By definition of the skew Schur functions and the Chern Simons invariants, 
SQ(q PM >v/„(<l Q+PM) ) = Nb^bq^' 1 ) = N Bh S b M Q' g ivin § the final result 



-Af|fi| 



Q- D,B,C 

where D is a representation of U(N), B and C are representations of U(M), and 

« = (-) |D| E<H H ^V ( 63 ) 
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Figure 4: The restricted region of general vertex describing U(N) x U (M) x U (K), for N = 3, 
M — 1, and K = 2, in the local C 3 geometry is demarcated by the green planes, and the 
asymptotes are labelled. 



5.3 The general vertex 

The most general possibility for intersecting toric D4 branes is the triple intersection 
with representations running along all three legs of the vertex. The partition function 
is significantly more complicated, and, although it can still be computed schematically 
with the transfer matrix technology, we will not be able to express it in terms of 
the Chern-Simons link invariants. In the context of OSV, the existence of a large 
N factorization remains clear, simply from the crystal picture itself, as explained in 
detail in the next section. We will also consider a special case that can be solved more 
explicitly. 

The partition function for the crystal oriented as in figure 4, but with M > N, can 
be determined using the same methods to be given by 



g -|ll^ll 2 -fllQII 2 g (^-^)(|fi|-|QI) x 



M-N 



< R\l[T ± ( q - n+1 2)v J] T±(q m - 1 *)ll r T(<i M ~ N+l ~ l2 )\Q t 



(64) 



n>0 m=l Z>0 



where we have chosen the framing with parallel rectangular walls for convenience. The 
pattern of signs of T is given by P, and the projection V forces the Young diagram to live 
in the restricted region determined by K in the vertical direction, and, in the horizontal 
direction, by the distance along the diagonal from the background representation, P, 
to the point (N,M). We orient the representations along the intersecting pairs of D4 
branes such that R is of type (N, K), Q is (K, M), and P is (M, N) in the notation of 

(25). 

Let us consider the special case that all of the boxes in the bifundamental region 
of P are filled, and joined to a U{N) representation, Pn, and a U(M) representation, 
Pmi as depicted in figure 5, for which a more detailed analysis is possible. This form 
of P implies that V = T 3 ^, and commuting the T±s appropriately, it is convenient to 
define 

W = q -h\\R t \?-m\\^-M)m-\Q\) X (P N )X{P M \ (65) 
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Figure 5: A special kind of (M, N) Young diagram is which the U(N) and U(M) degrees of 
freedom are decoupled is shown. 



in terms of the commutator term X (Pjv) =q 2 II- p atII 2 + 2 ' Pjv 'sp JV (g pN ), and the framing 
prefactor. Then we find that the statistical partition function is 

w E E s ^ f+w+ ^Kv^^^ 

L£U(K) \i,v 

(66) 

Now apply various identities to expand and simplify Schur functions evaluated at ex- 
pressions of the form q p+P N 5 to obtain, 

E E ^^ |i?h ^ lMl+(M+7V)|Lh ^ |!yh ^ IQI (-) IH+IH S/ ? / At (g pJV+P " ;v ) x 
LeUiK)w,a,p (67) 

s a /A^ +PN+PN ) s L/^ P ~ N )s L /p{q p - M )s 

Recognizing the appearance of the Chern-Simons link invariant from 1)53(1 . we can 
simplify the partition function to find 

W V E^ lAl+ ^ lQl (-) M+Wl+lAl+lBl s R/ll (qP N+P ~ N ) x 

A,B£U(K) w (68) 

Therefore putting everything together, and rewriting the result in terms of Chern- 
Simons invariants, we conclude that, 

Z = g (T- w )l«l+(T- JV )IO|-il|fl t ll 3 -^IIOII a g-ill^ll a -§l|Jl,ll a +^ \Pn\+%\Pm\ 

A,B<=U(K) 

C,DeU(N) E,FeU(M) 

(69) 

where the fusion coefficients are defined by equation (|63|). 
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6 Chiral factorization and conjugation involu- 
tion 



The chiral factorization of the indexed entropy at large D4 charge can be seen very 
elegantly in the melting crystals. The appearance of one factor of the topological string 
is immediate, since in the limit that N goes to infinity, the restrictions on the tableaux 
become vacuous, and we exactly reproduce the propagators and crystal vertices of the 
topological A-model in the same toric geometry. Referring to the derivation of the 
induced brane charges (|29|). we see that one obtains exactly the framing factors and 
vertex amplitudes of the topological A-model on the Calabi Yau [H]. Moreover, the 
extra contribution to the induced DO charge is precisely 

\[C]-[D], (70) 

the intersection number of the D4 branes on [C] with the homology class, [D], wrapped 
by the D2 branes. This combines with the chemical potential for the D2 charge to give 
the attractor value of the Kahler moduli, 

t a = \gto P [C] ■ [D a ] -i9 a , (71) 

for the curve [D a ] € H 2 (X). 

To see what the anti-chiral block arises from, recall that the propagator depends on 
the quadratic Casimir, so we expect both small representations and their conjugates 
to contribute, following 

C 2 {TZ) = k r+ + kr_ + N(\R+\ + \R_\) + M R 2 + 2£ R (\R+\ - |J2_|), (72) 

where 1Z = (RQ)[£r] is the U(N) representation with chiral part anti-chiral com- 
ponent and overall U(l) factor £r. 

There is an involution on restricted three dimensional partitions of the type we 
have been discussing which is inherited from the conjugation of U(N) representations. 
Consider two successive Young tableaux in the transfer matrix approach, R y Q, that 
is R\ > Qi > i?2 > Q2 > • • • • Then, when both are U(N) representations, under 
conjugation, one can see that Qn > Rn > . . . , and thus Q y R, up to a shift by a 
sufficiently large and positive U(l) factor. Therefore this defines an involution on the 
chiral fermion Hilbert space, sending 

\R >^ \R >, q Lo -» q~ Lo and T_ -» T+, (73) 

since T generates the interlacing relation, and |jF2| = — \R\. Naturally, these relations 
make sense only when we are restricted to U(N) representations. 

The involution also acts naturally when the slicing is oriented differently, for exam- 
ple, suppose R and Q are U(N) representations such that R l y Q l . It is easy to see 
that (R) l n + = N, for conjugation with shifting by some sufficiently large U(l) 

factor A. Thus (R) 1 -< (Q) , again up to an appropriate U(l) shift, chosen so that all 
column lengths are positive. 

Applying this transformation to the interlacing representations of the crystals, one 
can show that they are invariant. For a single stack of D4 branes, for example, we have 
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s{g = q ~\c 2{ R)-^ { Q) < m Yl r+iq-^vjj n r_(g™-S)|(Q)* > 

n>0 m>0 

= 9 -|^H^W)< jR *|J]r_(g-l)^ n r +(^ m+l )l^ >=5 S ) > 

n>0 m>0 

by commuting the T's in the obvious way, using crucially the fact that we are free to 
shift by ?7(1) factors to avoid having rows of negative length. 

Therefore the factorized form Zym ~ iV'top) 2 , follows automatically, by breaking all 
of the representations associated to D2 branes wrapping the legs into their chiral and 
anti-chiral pieces. The sum over chiral blocks (ie. ghost branes) first found in |3j arises 
because crystals with order N\A\ boxes can connect the chiral and anti-chiral regions 
with the ghost Young tableau A via a N x A transversal configuration, contributing 
e ~9sN\A\ ^ kg p ar titi on function, which remains finite even at large N. The simplest 
example of this is the vertex partition function of N coincident D4 branes, where in 
the large N limit diagrams, v, along the truncated direction, as shown in figure 1, 
are exactly the ghost tableaux, which clearly have no restrictions on their rank. The 
additional DO branes is the corner contributing to the chiral block will experience an 
effective background described by the topological vertex Cgt +R+U , as expected from the 
factorization formulae of and [I]. 



7 Concluding remarks and further directions 

We have solved the theory of BPS bound states of D2 and DO branes to D4 branes 
wrapping deformable cycles in toric Calabi-Yau by localization, obtaining a truncated 
crystal partition. This calculation is exact, even for small charges, and smoothly asymp- 
totes, in the OSV large charge limit, to the square of the topological string amplitude 
in the crystal vertex expansion. The crystal configurations we have found count the 
extra DO brane bound states at the singular points of the non-reduced surface obtained 
by giving nilpotent VEVs to the worldvolume adjoint fields controlling the normal de- 
formations of the D4 branes in the Calabi Yau. These are the toric invariant, zero size 
limit of the extra 2-cycles in H2(C) > H2(X), in the phase where the D4 charge results 
from a single brane wrapping the high degree, smooth surface C. 

The index of BPS bound states should be the same in all three cases, since it is 
protected by supersymmetry. We check that the results found here agree with the 
topological gauge theory calculations in known examples. This gives a nice realization 
of the underlying geometric nature of the interacting topological Yang Mills theories, 
since the crystal picture is not at all manifest in the quiver of four dimensional U(N{) 
gauge theories with bifundamental couplings along the intersections. 

We view this work as a step in the development of a nonperturbative topological 
vertex, that is a general method to compute the partition function of D4 branes on 
toric Calabi Yau manifolds in the OSV ensemble, which defines the nonperturbative 
completion of the topological A-model amplitudes. We have computed the relevant 
vertex and gluing rules, whose analogs in the topological string were sufficient to solve 
the theory, since the A-model only involves holomorphic maps from curves to the Calabi 
Yau. 
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The main problem in the exact OSV calculation on toric Calabi Yau which remains 
to be solved is the effect of branes wrapping compact 4-cycles. In toric Calabi Yau, com- 
pact 4-cycles always have positive curvature, and therefore negative self-intersection. 
As emphasized in [3], we must wrap D4 branes on very ample divisors to create black 
holes with large positive values for the geometrical moduli at the attractor point. This 
is the regime of validity of the OSV conjecture. Moreover, the theory living on D4 
branes instead wrapping rigid divisors with negative self-intersection cannot be solved 
by the Af = 1 mass deformation of [£], and it exhibits background moduli dependance 
and lines of marginal stability where the number of BPS states jump. 

This might seem to preclude the analysis of any more complicated geometries with 
multiple compact 4-cycles in the toric context, such as the resolution of A3 ALE fibred 
over P . However, if the D4 branes wrap toric 4-cycles, including some rigid cycles, 
which still sum to a very ample divisor, or equivalently, give rise to positive Kahler 
moduli at the attractor point, then there is no reason to expect a breakdown of the 
OSV conjecture. Therefore it seems likely that the effect of intersection with sufficiently 
many noncompact branes wrapping " good" 4-cycles will eliminate the problems of the 
twisted N = 4 theory on compact rigid surfaces. It would be very interesting to see 
how this works in detail, and to understand the new contributions arising from such 
branes. This could lead to a more complicated set of gluing rules, involving additional 
contributions from " loops" in the toric diagram as well as vertex and propagator terms 
that we have found here. 

It would also be interesting to investigate the relation between our crystals and 
those studied in |18j . which have a different kind of truncation, arising purely in the 
context of the q expansion of (perturbative) topological string theory. This might 
further elucidate why the Chern Simons amplitudes S^j appear, even at finite N, in 
the solution of D4 brane U (N) gauge theories. 

Acknowledgements 

I would like to thank Mina Agaganic, Joe Marsano, Natalia Saulina, Xi Yin, and 
especially Cumrun Vafa for valuable discussions. This work was supported in part 
by NSF grants PHY-0244821 and DMS-0244464. I would like to thank the Simons 
Workshop at the Yang Institute of Theoretical Physics at Stonybrook, where part of 
this work was conducted, for hospitality. 

Appendix A 

All r + commute, as do T_, however one can show that 

r_(x)r + (y) = r + ( y )r_(x)(i-^). (75) 

Also, 

T±(x)q nL ° = q nLo T±(xq- n ). (76) 
These operators are conjugate; following the definition (|4Ujl it is easy to see that 

(r + (x))t = r_(x- 1 ). (77) 
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(80) 



The repeated application of the operators results in the appearance of skew Schur 
functions, since 

Y[T-{x n )\u>= J2 s n/v( x n)\v >■ (78) 

We will often come across skew Schur functions evaluated at special values of the form 
%n = Pn + v n , where p = {— ^, — |, . . . }, and v is a Young diagram, as well as the 
finite N version, X{ = pf + R{, i = 1, . . . N, where R is a U(N) representation, and 
pf = — ^ is the Weyl vector. These functions obey many beautiful relations, 
including 

where Nj^ are the tensor product coefficients. There are various useful summation 
formulae as well, 

S n/ri( x ) s v/r)(.y) = Yl( 1 ~ x nVm) ^ S x / ^{v) s \ / v{ x ) 
r] n,m X 

^2s^ /u {x)s u {y) = s u {x,y), 

V 

where (x, y) denotes the conjunction of the strings of variables x and y. 
We also make use of the quadratic Casimirs, 

Kfi = H/ijl 2 - H^ll 2 , and C^ N) (R) = N\R\ + k r , (81) 

for any Young diagram \x and U(N) representation R. 
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